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Extending work of E, Khanin, Mazel and Sinai jl| on the one-dimensional 
Burgers equation, we show that density pdf's have universal power-law 
tails with exponent -7/2. This behavior stems from singularities, other 
than shocks, whose nature is quite different in one and several dimensions. 
We briefly discuss the possibility of detecting singularities of Navier-Stokes 
turbulence using pdf tails. 

1 Introduction 

In recent years there has been considerable interest in probability density func- 
tions (pdf) for Navier-Stokes turbulence. Similar questions can be asked for 
random solutions of Burgers equation ("burgulence"). We are interested here in 
the tail behavior of the pdf of the density p for solutions to the d-dimcnsional 
Burgers equation in the limit of vanishing viscosity (y — > 0): 



The initial potential ^o(ro) and the initial density po(ro) are random functions of 
the space variable. This problem arises, for example, in the study of large-scale 
structures in the Universe (see Ref. Q and references therein). As is well known, 
the Burgers equation leads to shocks in which the density of matter is infinite. 
Yet, large but finite densities do not necessarily occur in the neighborhood of 
shocks. E et al. considered a related problem of determining the pdf of 
the velocity gradient for the one-dimensional Burgers equation with large-scale 
and white-in-time random forcing. They showed that large negative gradients 
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Figure 1: Pdf of the velocity gradient at negative values in log-log coordinates, 
for one-dimensional kicked burgulence. A power law with exponent —7/2 is 
obtained over two decades. Simulation uses a modification of the Fast Legendre 
Transform algorithm of Ref. (|] . 



come from "preshocks" (nascent shocks) which contribute a power-law tail with 
exponent -7/2 to the pdf. Preshocks correspond to fast fluid particles catching 
up for the first time with slow ones. They constitute discrete events in (Eulcrian) 
space time. We have shown that the -7/2 law for the pdf of (negative) velocity 
gradients in one dimension applies also (i) for decaying (unforced) "burgulence" 
with smooth (i.e. large-scale) random initial conditions [g| and (ii) for the case of 
a deterministic time and space-periodic force which is a sum of delta functions 
in time with smooth space dependence [Q. For this "kicked burgulence", the 
-7/2 law was also obtained numerically (see Fig. 1). 

2 "Kurtoparabolic" points and the -7/2 law 

We turn now to the multi-dimensional problem (|l|)-(^|) with smooth and ran- 
dom initial conditions (for details see Ref. Q). As is known, the one or multi- 
dimensional viscous Burgers equation can be solved explicitly by means of the 
Cole-Hopf transformation. From this, by taking the limit v — > 0, a "maximum 
representation" can be derived for the velocity potential 

</>(r, t) = max (VoM - ^T^) ■ (3) 

The maximum is achieved at a point at which the gradient of the r.h.s. vanishes, 
leading to r = rg + £v(ro,0). From this it follows that Tq is a Lagrangian 
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coordinate. The "naive" Lagrangian map ro i— > r given by this relation is however 
not invertible, except for short times. The requirement that the maximum in 
(Q) is global can be recast in geometrical terms by introducing the Lagrangian 
potential <p(tq, t) = tip (r ) — |r | 2 /2 and the "proper" Lagrangian map r i— > r = 
— Vy> c (ro,t), where tp c is the convex-hull of ip with respect to ro. For example, 
in one dimension, the graph of ip c is obtained by tightly pulling a string over 
the graph of ip. The graph of (p c coincides with that of if at regular points, 
wherever fluid particles have not yet fallen into shocks. It also contains linear 
and ruled manifolds associated to the different types of shocks: segments when 
d = 1, triangles and ruled surfaces when d = 2, etc. Conservation of mass 
implies that the density is given, at regular points, by p{r,t) = po( r o)/^( r o, t), 
where J is the Jacobian of the Lagrangian map. (The density is infinite in 
shocks.) Since the Jacobian is (up to a factor (— l) d ) equal to the Hessian of the 
Lagrangian potential (determinant of the matrix of second space derivatives), 
it follows that large densities are typically obtained only near parabolic points. 
However, arbitrarily close to a parabolic point there are generically hyperbolic 
points where the surface defined by ip crosses its tangent (hypcr)plane and which, 
therefore, do not belong to its convex hull. Yet, there exist in general exceptional 
"kurtoparabolic" points which are parabolic and belong to the boundary of the 
set of regular points. Near such points, arbitrarily large densities are obtained. 
In one dimension, the only kurtoparabolic points are the preshocks which are 
discrete space-time events in both Eulerian and Lagrangian coordinates. In 
two and more dimensions, kurtoparabolic points are also born at preshocks but 
persist in general for a finite time (see Fig. 2). In Eulerian space, they are 
associated to boundaries of shocks (e.g. end points of shock lines for d = 2). 




Figure 2: The Lagrangian potential in two dimensions, just after a preshock 
(left) and in the neighborhood of a kurtoparabolic point (right). Continuous 
lines: separatrices between the regular part and the ruled surface of the convex 
hull; dotted lines: vanishing of the Jacobian of the Lagrangian map. A and A' 
are a pair of kurtoparabolic points born with the shock. 

We just indicate general ideas involved in the determination of the density 
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pdf for random initial conditions (which need not be homogeneous). The de- 
termination of the large-p tail of the cumulative probability distribution of the 
density is equivalent to finding the fraction of Eulcrian space-time where p ex- 
ceeds a given value. The latter is determined by changing from Eulerian to 
Lagrangian coordinates and using a suitable normal form (i.e. a Taylor expan- 
sion to the relevant order) of the Lagrangian potential near a kurtoparabolic 
point. The theory is rather different in one dimension and higher dimensions, 
because kurtoparabolic points are persistent only in the latter case. However, 
the scaling law for the resulting pdf, namely cx p~ 7 ^ 2 , is the same in all dimen- 
sions. In fact, when d > 2, two orthogonal spatial directions play the same role 
as space and time in one dimension. 



3 Detecting Navier— Stokes singularities 

It is clear that the algebraic tail of the pdf of velocity gradients or of density for 
burgulencc comes from identified singularities. Measurements of pdf 's for space 
or time derivatives of Eulcrian velocities for incompressible three-dimensional 
Navier-Stokcs turbulence have not revealed power law tails, but such tails may 
just have been, so far, "lost in the experimental noise" . There has indeed been 
considerable speculations about singularities of the Navier-Stokes equations in 
the inviscid limit. If singularities with divergent gradients are present, they will 
give power-law tails, at least as intermediate asymptotics when the viscosity is 
small (the converse is however not true, since statistical effects not related to 
singularities can also give power laws). The confirmed absence of power laws 
would probably rule out singularities. 
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